AD-A203  698 


I  inC  FILE  w 


INVESTIGATION  OF  EXTENSIONS  TO  THE  DISTORTED  BORN 
APPROXIMATION  IN  STRONG  FLUCTUATION  THEORY 
FINAL  REPORT 


CONTRACT  NO.  NOOOW-87-0784 


Prepared  For 

OFFICE  OF  NAVAL  RESEARCH 
DEPARTMENT  OF  THE  NAVY 
Arlington^  Virginia  22217 


Report  No.  9316 


STATEAfENT  A 
Approved  for  public  reload 
L  Djgtribution  Unlimited 


DTIC 

ELECTE 
JAN  3  1)989 


October  1988 


INVESTIGATION  OF  EXTENSIONS  TO  THE  DISTORTED  BORN 


APPROXIMATION  IN  STRONG  FLUCTUATION  THEORY 

Final  Report 

October  1988 
Report  No.  9316 

-C- 

Contract  No.  N00014-87-0784 

by 

A.  Stogryn 

Prepared  for 

OFFICE  OF  NAVAL  RESEARCH 
DEPARTMENT  OP  THE  NAVY 
Arlington,  Virginia  22217 


Prepared  by 


AEROJET  ELECTROSYSTEMS  COMPANY 
1100  West  Rollyvale  Street 
Azusa,  California  91702 


Report  9316 


INTRODUCTION  AND  SUMMARY 

This  study,'  funded  under  Contract  N00014-87-0784 considers  some 
scattering  problems  relevant  to  the  analysis  of  brightness  temperatures  and 
scattering  cross  sections  of  sea  ice  and  its  snow  cover.  It  is  a  continu- 
ation  of  work  performed  at  Aerojet  ElectroSystems  Company  under  previous 
contracts  with  the  Office  of  Naval  Research.  The  theme  of  these  studies 
has  been  the  application  of  strong  fluctuation  theory  to  the  description  of 
the  microwave  radiometric  properties  of  random  media  such  as  sea  ice  and 
snow. 


These  past  studies  have  demonstrated  the  applicability  of  the  bilocal 
approximation  to  the  equations  of  strong  fluctuation  theory  in  calculating 
electromagnetic  material  properties  of  random  media  such  as  effective  di¬ 
electric  constants  [1],  (2J.  The  resulting  equations  have  received  con¬ 
firmation  from  numerous  experiments  (see  (31,  (4)  for  measurements  on  snow 
and  the  discussion  in  [2]  for  experimental  data  on  sea  ice).  To  the  pre¬ 
sent  time,  these  are  the  only  published  equations  accounting  for  the  order 
of  magnitude  increase  in  the  dielectric  loss  experienced  by  snow  at  fre¬ 
quencies  above  20  GHz  as  compared  to  its  low  frequency  value. 

The  bilocal  approximation  has  also  been  used  in  conjunction  with  the 
distorted  Born  approximation  in  calculating  the  brightness  temperature  of 
snow  [3]  and  sea  ice  [5],  The  distorted  Born  approximation,  which  uses 
only  the  mean  electric  field  within  the  random  medium  to  evaluate  scat¬ 
tering  and  thermal  emission  into  the  air  above,  has  had  some  success  in 
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explaining  measurements  of  the  brightness  temperature  of  snow  at  fre¬ 
quencies  below  37  GHz,  young  and  first  year  ice  at  frequencies  up  to  at 
least  140  GHz,  and  old  sea  ice  at  frequencies  into  the  20  GHz  range.  Its 
use  in  conjunction  with  strong  fluctuation  theory  has  led  to  results  con¬ 
siderably  superior  to  that  found  by  the  application  of  radiative  transfer 
theory.  However,  for  media  such  as  snow  and  old  sea  ice,  the  distorted 
Born  approximation  leads  to  predictions  of  emissivities  and  brightness 
temperature  which  are  considerably  higher  than  observed  values  at  suffi¬ 
ciently  high  frequencies.  The  reason  for  this  behavior  is  suspected  to  lie 
in  the  neglect  of  scattering  by  the  random,  or  incoherent,  field  within  the 
medium  when  calculating  the  emissivity  and  scattering  coefficients. 

■  a  / 

It  is  the  purpose  of  *the  work  ^reported  here  to  examine^how  -thi^  here- 
tofore  neglected  contribution  of  the  incoherent  field  to  the  emissivity  and 
scattering  cross  sections  of  a  random  medium  may  be  taken  into  account. 

^Ther«-axe-  two  main  technical  problems  are  solved.  In  Section  II, 

strong  fluctuation  theory  equations  are  derived  for  the  second  moments  of 
the  electric  field  in  an  anisotropic  medium.  Previous  work  in  this  area 
had  not  considered  anisotropy  in  the  random  medium  so  that  it  was  not  clear 
what  form  the  second  moment  equations  would  take  in  the  case  of  a  general 
random  medium  such  as  sea  ice,  where  anisotropy  may  occur.  Nor  was  the 
most  useful  explicit  equation  written  for  the  isotropic  case.  The  analysis 
is  dependent  on  obtaining  equations  which  allow  the  elimination  of  the  aux¬ 
iliary  field  jg,  which  always  arises  in  strong  fluctuation  theory,  in  favcr 
of  the  electric  field  g,  which  is  easy  to  work  with  when  electromagnetic 
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boundary  conditions  are  applied  and  which  has  a  more  direct  physical 
meaning  than  F. 

Section  III  of  this  report  makes  use  of  a  . particular  equation '^derived 
in  Section  II  as^he  basis  for  a  detailed  reduction  of  the  second  moment 
equation  to  a  form  which  is  practical  for  numerical  computation.  The 
problem  of  a  layered  medium  with  a  spherically  symmatric  correlation  func¬ 
tion  is  considered.  The  spherical  symmetry  of  the  correlation  function 
implies  that  a  scalar  description  of  the  dielectric  properties  of  the 
medium  is  valid  and  results  in  the  elimination  of  several  technical 
problems  which  would  arise  had  an  anisotropic  correlation  function  been 
assumed.  Thus,  the  theory  is  applicable  to  snow  cover  on  sea  ice  as  well 
as  old  sea  ice  (provided  that  the  correlation  functions  describing  brine 
pockets  in  the  ice  is  spherically  symmetric  as  is  suggested  by  some 
studies)  but  will  not  be  directly  applicable  to  young  and  first  year  ice 
which  are  anisotropic.  However,  the  development  of  the  theory  is  general 

- A  ) 

enough  to  allow  complex  vertical  variations  in  the  mean  dielectric  profile 
of  the  medium  so  that  it  may  be  applied  to  realistic  layered  structures 
found  in  naturally  occurring  snow  and  ice.  Particular  care  is  taken  to 
identify  and  treat  all  singular  behavior  which  could  cause  problems  in  a 
numerical  solution. 

The  equations  developed  in  Section  III  have  not  yet  been  programmed 
for  solution  on  a  digital  computer.  In  view  of  the  expectation,  based  on 
good  physical  arguments,  that  the  use  of  these  equations  will  greatly 
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extend  the  high  frequency  range  of  validity  of  strong  fluctuation  theory 
over  the  distorted  Born  approximation  for  strongly  scattering  media,  it  is 
recommended  that  the  next  phase  of  these  studies  be  devoted  to  writing  a 
FORTRAN  program  to  solve  these  equations.  There  also  exists  the  possi¬ 
bility  that  one  of  the  integrations  involved  in  solving  these  equations  may 
be  performed  analytically,  thus  substantially  decreasing  the  complexity  of 
the  proposed  computer  program.  Hence,  an  additional  analytic  effort  in 
this  area  would  be  well  worthwhile.  It  is  further  recommended  that  the 
computer  program  that  is  written  be  applied  to  examining  the  brightness 
temperature  behavior  of  an  isothermal,  dry  snow  pack  at  microwave  fre¬ 
quencies  where  the  distorted  Born  approximation  begins  to  fail.  Compa¬ 
risons  with  experimental  data  should  be  made.  There  is  sufficient  experi¬ 
mental  data  available  to  make  such  comparisons  enlightening.  It  should 
also  be  useful  to  see  what  these  new  equations  predict  for  the  brightness 
temperature  of  old  sea  ice  (at  a  uniform  temperature)  assuming  that  the 
applicable  brine  pocket  correlation  function  is  spherically  symmetric. 

The  effect  of  snow  cover,  which  usually  occurs  on  the  ice,  should  also  be 
studied. 
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STRONG  FLUCTUATION  THEORY  EQUATIONS  FOR  ELECTRIC 
FIELD  SECOND  MOMENTS  IN  ANISOTROPIC  MEDIA 


I .  INTRODUCTION 

Strong  fluctuation  theory  for  random  media  ([ll-[5])  has  found 
application  to  numerous  interesting  geophysical  problems.  In  [6],  it  was 
shown  how  the  equations  for  the  mean  electric  field  in  the  bilocal  approxi¬ 
mation  and  the  random  field  in  the  distorted  Born  approximation  could  be 
formulated  directly  in  terms  of  the  electric  field  instead  of  in  terms  of 
equations  containing  an  auxiliary  field  F  (see  [l]-[4],  [6]).  BreJchovskiJch 
(7]  has  shown  how  equations  for  electric  field  second  moments  could  also  be 
obtained  in  strong  fluctuation  theory  in  the  case  where  the  fluctuations  in 
the  random  medium  have  a  spherically  symmetric  correlation  function  and 
hence  lead  to  a  scalar  rather  than  a  second  rank  tensor  description  of  the 
quantity  ^(  [ 1 ]-[ 6 ] )  which  is  used  to  describe  the  fluctuations  of  the  medium 
in  strong  fluctuation  theory.  In  this  communication,  the  equations  for  the 
second  moments  are  generalized  to  the  case  of  an  anisotropic  medium. 

II.  SECOND  MOMENT  EQUATIONS 

The  electric  field  E(r)  in  a  random  medium  described  by  a  random 
second  rank  dielectric  tensor  K  is  governed  by  the  equation 

(L  +AL)  E  =  0  (1) 

o  ~ 

where  L^+AL  is  the  operator  VX  VX  -  •'p  is  t*'®  free  space 
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propagation  constant.  In  strong  fluctuation  theory  a  non-random  quasi¬ 
static  dielectric  tensor  K  is  introduced  in  terms  of  which  the  operators 

o 

and  A L  are  given  by 


(L 


)  .  - 
O  2.J 


8  X.  8x  . 
1  3 


6.  .V^- 

13 


O  O  1] 


(2) 


AL=  k  ^(K  -K*^) 
o  o 


(3) 


The  quasi-static  dielectric  tensor  is  'defined  by  means  of  the  equation 


<  4  >  =  0 


where  the  angular  brackets  denote  expected  value  and 


^  =  AL[1+S^] 


(4) 


(5) 


Here,  all  second  rank  tensors  are  represented  by  matrices  and  multipli¬ 
cation  indicates  matrix  multiplication.  The  matrix  S  in  (5)  is  defined  as 
the  coefficient  of  the  delta  function  in  the  Green's  function  T  determined 
by  (2)  which  may  be  expressed  as 


r(r,  r')  «  S6  (r-r’)  +  pvr’(r,r') 


(6) 


In  this  equation,  PV  means  that  the  principal  value  is  to  be  taken  in 
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integrals  in  which  T'  appears,  with  these  definitions,  an  auxiliary  field 

P  =  [l+SALl  E  (7, 

is  introduced  in  strong  fluctuation  theory.  It  has  the  property  that  it 
satisfies  an  integral  equation  whose  kernel  is  free  of  delta  functions, 
namely 

F  =  G  4  F  (8) 

where  G  is  the  operator  whose  kernel  is  PVF’ .  Equations(5)  and  (7)  show 
that 

^  F  =  Ale  (9) 

•  ^  ^ 

while  (7)  and  (9)  show  that 


F  =«  E  +  S  £  F 


(10) 


or 

E  =  [1-  S  4  ]  F 


(11) 


In  the  following  derivation  of  the  second  moment  equations,  the 
argiunents  of  the  various  quantities  involved  will  be  abbreviated  to 
subscripts.  Thus,  e.g.,  the  notation 

'ill  S‘<£2»  (121 


is  introduced  where*  indicates  complex  conjugate  and  the  juxtaposition 
of  two  vectors  is  simply  the  matrix  representing  the  direct  product  of  the 
vectors.  Following  (7],  it  is  convenient  to  define  the  operators  ^  and  B 
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by  the  equations 


4,  <F,  F  *>  =  <  4,P,  F  *: 


:i3) 


^  <£l  l2^  =  <£l  ^2*  £2*" 


(14) 


C,  <F,  F  *>  =  <F,  S  •  £  *  F  •> 


(15) 


B  <F,  F  *>  =  <  £,  F,  £  *  F  *> 


(16) 


B,<  F,  F  *> 
1  ^2 


=  <  ^1  ^1  £1^2*  r 


(17) 


B,<F,  F  *  >  =  <  £  ,  F,  S  *  I  *  F  *  > 

2  ^2  1  2  2  •n»2 


(18) 


Equations  (15),  (17),  and  (18)  are  additional  definitions  that  are  required 
compared  with  [7]  because  S  does  not  commute  with  the  fields  in  general. 

For  the  special  case  of  a  spherically  symmetric  correlation  function  for  | 
and  scalar  (and  hence  scalar  /  S  is  simply  the  scalar  [8J 


S  -  -1/(3  k  K  ) 
o  o 


(19) 


and  may  be  removed  from  the  angular  bracicets  in  these  equations. 


The  operators  ^  and  B  allow  a  relation  between  <F,  F  * >  and  <E,  E  *: 

«*1  <^2 
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to  be  established.  A  calculation  using  (7),  (9)  and  (11)  yields 


h  £2*=  Si  S2*  f^'^l  tl^  Si  ^2*  ^2*  I2*  ^  ®1  ^1  Si  ^2*’S2* 


^  ^  Si  ^2*  §2*  S2* 


(20) 


which  has  the  consequence 


-1 


<Jj  £/>  -  [1-  43-  Sj  ♦  Sj  <£3  E3*> 


(21) 


Now  ( 9 )  and  ( 1 )  show  that 


£  F  fi  *  F  *»ALE  AL  *  E  * 

^1  Si  62  S2  ~1  S2 


=  L  E  L  *  E  * 
ol  Si  o2  S2 


(22) 


so  that,  using  (16)  and (21),  the  equation 


‘■01  K2‘  - 


-  "1^1  *  "i  “i'M'Si  £2* 


>  -  0 


(23) 


is  found.  Here  L  _*  is  defined  as 
02 


^^o2*^ij)ci  "^Sl  £2  *^1)1  ^^o2  \  l  "^Sl  £2 


(24) 
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where  the  sununation  convention  is  used  for  repeated  subscripts.  Another 
interesting  equation  is  obtained  from  the  pair  of  equations 


h 


(25) 


•,  F  =  AL,  E,  F.,*=  AL,  E,  [E. 

1  >*.1  ~2  1  '~2  i  "•!  ~2 


+  S, 


£2 


(26) 


which  follows  from  (17)  and  (9).  Subtracting  (25)  from  (26),  using  (1),  and 
taicing  expected  values  yields 


L,  <E  E  *>  +  [4,  -  B  ]  <F  F  *> 
ol  ~1  ~2  12  ~1 


(27) 


Of  course,  the  substitution  of  (21)  in  (27)  yields  an  equation  containing 
<E,  E  * >  alone.  The  equations  (again  following  from  (7)  and  (9)) 


h  \  h  ^2*  !:2‘ 


Si  A 


£1 


6L. 


E,' 

'~2 


(28) 


F,  i*  F  *  »  F,  AL  *  E  ‘ 
^1  ^2  ^2  2  ^2 


[E,  +  S,  A  L,  E,  ]  AL  *  E,* 
~1  1  1  ~1  2  0-2 


(29) 


yield  an  equation  analogous  to  (27): 


<E-  E  *>  +  [  4,-  Bj  <F  F  *> 
«•!  '^2  2  1  *^1  **•2 


0 


(30) 


III.  APPROXIMATIONS 

Equations  (23),  (27),  and  (30)  are  exact.  If,  as  in  [6],  quadratic 
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and  higher  order  terms  inland  B  are  ignored,  then  the  equations 

Kl  ‘o2‘-  ”^£1  52*’=  “ 

j‘'ol  "  “2f  '‘2*^  -o 

)l  5*+  c,-  B  i  <E  B^*>=  0 
(  02  2  1)  ~i 

are  found.  Allowing  for  slight  differences  in  notation,  eq  (32)  reduces  to 
that  given  in  [7]  (written  for  a  vanishing  external  current  density)  if  the 
medium  is  isotropic  so  that  (19)  may  be  used. 

A  slightly  different  version  of  (31)  is  useful  in  problems  where  the 
mean  field  <E>  has  already  been  calculated  and  an  expression  for  the  inco¬ 
herent  contribution  to  <E,  E.*>  is  needed.  That  is,  if  the  field  is 

•*■1 

written  as 


(31) 


(32) 


(33) 


E  ■»  <E>  +  E*^ 
^  ^ 


(34) 


where  E^  is  the  random  part  of  E,  then  <E,^  E_^*>  is  to  be  found.  Perhaps 
the  easiest  way  to  derive  the  desired  equation  is  to  note  that,  according 
to  [6],  E^  satisfies 


-  <ep>f 


(35) 
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exactly .  Then 


'£/  £2"*^-  ^02''  i“  <  £1  £2'^-  <h  £1^  ‘«2*  £2*4 


(36) 


Now,  for  large  separations  r -r.,  (16)  shows  that  B  <F  F  *> _ , 

<|j^  F^>  <^2*  ~2*^  that  the  bracketed  quantity  in  (36)  vanishes,  on  the 

other  hand,  for  small  to  moderate  values  of  r  -r_,  the  first  term  in  the 

'“I  ~2 

bracket  dominates.  Thus,  if  it  is  assumed  that  the  field  at  a  point  is 
only  slightly  dependent  on  the  fluctuation  at  that  point,  (16)  yields 


<£i 


=  <‘=12’i,lJV2’> 


ki 


(37) 


where 


(C'  )  .  ..  ,  =  <  (C,  )  u  •*> 

12  ijki  ^1  ik  ^2  3I 


(38) 


is  the  4th  rank  covariance  tensor  for  the  random  tensors  ^  and  Since 

(38)  vanishes  for  large  r  -r.,  a  reasonable  approximation  to  the  bracket  in 
(36)  is  simply  (37).  Further,  since  higher  order  terms  are  already 
ignored  in  (31),  <P,  F_*>  in  (37)  may  be  replaced  by  <E,  E_*>.  Thus 


<E,^  E-^*>  -  L  C'  <E,  E,*> 

M.1  ^*2  ol  o2  12  'tel 


(39) 


Upon  noting  that  (34)  yields 


<Ej^  *  <  E  <  E  >2  + 


(40) 
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(35)  may  be  replaced  by 


<E>,  <E*>^ 
~  1  "  2 


(41) 


Either  (39)  or  (41)  represent  a  suitable  formulation  for  practical  compu¬ 
tation. 
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SECTION  2 

STRONG  FLUCTUATION  THEORY  CALCULATION 
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STRONG  FLUCTUATION  THEORY  CALCULATION  OF  SCATTERED 
FIELDS  FROM  MULTI-LAYERED  RANDOM  MEDIA 


I .  INTRODUCTION 


The  theory  of  electromagnetic  scattering  and  microwave  thermal 
emission  from  random  media  such  as  snow,  sea  ice,  and  vegetation  has  become 
an  important  area  of  research.  One  approach  to  studying  this  problem  has 
been  strong  fluctuation  theory  (1],  [2]  which  has  been  applied  to  calcu¬ 
lating  effective  dielectric  constants  [3]  of  snow  and  sea  ice  [4] -[6]  (with 
strong  experimental  support  [7] ) ,  scattering  cross  sections  of  vegetation 
[8],  [9]  and  snow  flO],  and  the  microwave  brightness  temperature  of  dry  and 
wet  snow  {3].  In  applications  to  scattering  and  brightness  temperature 
calculations,  the  distorted  Born  approximation  has  been  used  in  the  past. 
This  approximation  neglects  scattering  of  the  incoherent  field  and  results 
in  an  underestimate  of  the  true  scattering  by  the  medium.  Its  effect  is 
particularly  noticeable  at  higher  frequencies. 

It  is  the  purpose  of  this  work  to  improve  upon  the  distorted  Born 
approximation  in  strong  fluctuation  theory.  Integral  equations  are 
developed  for  the  second  moments  of  the  incoherent  electric  field  in  the 
case  of  a  multi-layered  random  medium,  each  of  whose  layers  is  bounded  by 
plane  surfaces.  It  is  shown  how  the  solution  of  these  equations  is  incor¬ 
porated  into  the  expression  for  the  bistatic  scattering  cross  section  of 
the  medium.  Of  course,  by  use  of  the  results  of  Peake  [11],  these 
scattering  cross  sections  may  be  used  to  calculate  the  thermal  emissivity 
of  the  medium  if  it  is  at  a  uniform  temperature. 
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II.  BISTATIC  SCATTERING  COEFFICIENTS 

Consider  a  layered  medium,  as  is  illustrated  in  Fig.  1,  that  has  a 
random  scalar  dielectric  constant  K^(r)  dependent  on  the  point  r  for  z<0. 

It  is  assumed  that  the  expected  value  of  K  ^ is  a  function  of  z  only  and  that 
the  covariance  function  of  is  spherically  symmetric  (i.e.  a  function  of 
the  magnitude  of  the  separation  of  two  points  only).  Further,  assume  that 
changes  in  the  expected  value  of  over  a  correlation  length  are  small 
(except,  of  course,  at  possible  boundaries  where  discontinuities  of  the 
dielectric  properties  arise). 

Suppose  that  the  medium  is  illuminated  by  an  incident  plane  wave 

E  =  E  n  exp(ik  •  r]  (1) 

««in  o  £o  ~ 

where 

Jt  =  )c  (sin0  ,O,-cos0  )  (2) 

~o  o  o  o 


and  p  is  a  unit  polarization  vector  which  will  be  taJten  to  be  either  hori- 
zontal  (h)  or  vertical  (v); 


A 


A 


h 


or 


A 

V 

— o 


10, -1,0) 


(cos  9  ,0,sin9  ) 
o  o 


(3) 


2-3 


Report  9316 


We  are  interested  in  the  scattered  power  contained  in  a  small  solid  single 
about  an  arbitrary  propagation  vector 

k  *  k^(  sin  9  cos4), sin0sin<i),cos 0)  (4) 

The  electric  field  which  is  generated  by  the  wave  (1)  may  be  expressed 
as 

E  =  £”*+  (5) 

where  E^'is  the  mean  wave  and  E^  the  random,  or  incoherent  part.  Using  the 
bilocal  approximation  of  strong  fluctuation  theory,  E™  satisfies  the  equa¬ 
tion  112] 

|l^-<5GC>  j>  e"'=  0  (6) 

where  is  the  operator 

Lq  -  37x7x  -  kg^  if  z>0  (7) 

Wx7x  -  k  if  z<o 

o  o 

and  G  is  the  non-delta  function  part  of  the  dyadic  Green's  function  £  (r,r' ) 
for  (7).  Here,  angular  brackets  are  used  to  denote  expected  values  and  is 
the  quasistatic  dielectric  constant  of  the  medium  and  is  defined  by  the 
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equation 


<C>  =  0 


(8) 


where 

5  =  AL[1+SAL]  ^  (9) 

and 

AL  =  k  ^  (K  -K^)  (10) 

o  o 


2 

In  (9),  S»  -l/(3k  K  )  is  the  coefficient  of  the  delta  function  part  of  ^  . 

O  O  s: 

The  second  moment  dyad  of  the  incoherent  field  satisfies  the  equation  [13] 

<E'^(r  )E^*(r  )  >=  J  d^r 'd^r  "C'(  1  ]  )£  )<E(r' )E*{r"  >>r(r  ,r"  ) 

^12  ^ 

V. 

where 

c'(|r'-r"J)«  <^(r' )  5*(r")>  (12) 


is  spherically  symmetric  according  to  the  assumptions  made  above.  Here  * 
indicates  complex  conjugate,  and ^ indicates  adjoint.  Notice  that  (5)  implies 
that 


When  this  is  used  in  the  right  hand  side  of  (11),  it  is  seen  that  (11)  is  an 
integral  equation  for  <E*^(  r^^  )E^*(  r2)> . 


The  bistatic  scattering  coefficients  defined  in  [11]  require  the 
decomposition  of  the  left  hand  side  of  (11),  evaluated  at  r^*  r^*  r,  into  an 
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angular  spectrum.  Namely,  the  components  of  (11)  are  written  as 


<E^(r)Ej^*(r)>= Jdfi<E^(k)E  r*(k)> 


where  d is  an  element  of  solid  angle  about  the  vector  k  and  the  integral  is 
taken  over  the  upper  hemisphere.  In  (14),  E  (k)  denotes  the  Fourier 
transform  of  E  (r).  The  scattering  coefficients  are 


Y  .  (k  k)=  (4Trcos9/E  ^cosO  )b.b.<E^(k)E.^*(k)>  (15) 

ab  ~o/~  o  oi]i"'3~' 

A  A 

where  a  and  b  denote  unit  vectors  describing  the  polarization  of  the  incident 
and  scattered  waves  respectively.  We  will  be  concerned  with  the  case  where 

A 

is  given  by  (3)  and  b  represents  horizontal  or  vertical  polarization  for  the 
propagation  vector  ( 4 ) s 

A 

h»  ( sin(fr, -cos  <^,0) 

or  ( 16 ) 

( -cos0cosip,-cos9sin<t>,sin0) 

The  computation  of  E*”  for  a  layered  medium  of  the  type  oeing  con¬ 
sidered  here  has  been  considered  in  detail  previously  [5],  [14].  These  works 
also  provide  detailed  expressions  for  <E^(k)E,*^*  (jS)>  in  the  distorted  Born 
approximation  which  is  obtained  by  neglecting  the  term  <E^( r^ )E^* ( r" ) >  when 
(13)  is  substituted  into  the  right  hand  side  of  (11).  The  task  here  is  to 
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retain  this  term  in  (11).  Thus,  write  (15)  as 


B  I 

u  (K.  k)-hY  .Ik  k) 

ab  — o,~  ab  ~o,—  ab  ~o,~ 


(17) 


B 

where  Y  represents  the  distorted  Born  approximation  contribution  to  Y  and 
do  do 

may  be  considered  to  be  )tnown.  The  contribution  of  the  incoherent  field  Y 

ab 

arises  from  the  second  term  on  the  right  hand  side  of  (13)  when  substituted 
into  the  right  hand  side  of  (11)  and  yields,  using  the  notation  for  the  two 
dimensional  Fourier  transform  of  the  Green's  function  developed  in  [14], 


Y  ()c  »  )c)=»  [)c  cos9/(7rcos  9  )]b,  b.  j  dz' 
ab~o  -  o  o  1  j 


dz"  A  (  z'  ,k  )c  ) 

i)c  X,  y 


(18) 


A*  (z"  ,k  ,k  )P^„(k  k  z',z"  ) 
J£  X  y  kZ  y,  y, 


where  the  dyad  F  is  the  two  dimensional  Fourier  transform 


)exp{i)t  x+i)c  y] <E*^(  r '  )E  ^ (  r"  )  >  (19) 

y  0^  ^0 

and  the  summation  convention  is  used  for  repeated  Latin  letter  subscripts. 

In  (19),  in  addition  to  the  z'  ,z''  dependence,  the  expected  value  on  the  right 
hand  side  is  a  function  of  x=  x' -x'  and  y=  y  -y"  because  of  the  horizontal 
translational  invariance  of  the  problem  under  consideration. 


F()c  k  z',z")=  dxdy  c'(|r'-r" 
s  y,  y,  ^  I 


Considerable  simplification  in  further  calculations  may  be  achieved  if 

use  is  made  of  an  additional  symmetry  of  the  problem.  Notice  that  the 

m  m* 

driving  term  in  (11)  is  E  E  which,  in  turn,  is  determined  by  the  incident 

wave  term  E,  (r  )E.*(r^)  arising  from  (2).  But,  for  either  horizontal  or 
~in  ~1  ~in  ~2 

vertical  polarization,  this  dyad  is  invariant  under  reflections  in  the  x,  z 
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plane.  The  geometry  illustrated  in  Fig.  1  as  well  as  the  expected  value  of 
the  fluctuations  is  also  invariant  under  this  reflection.  Hence,  so  is 
<e'^(  r,  ( r^  )  >  .  But  this  implies  that  the  (x,y),  (y,x),  (y,z)  and  (z,y) 
components  vanish.  Hence,  we  may  write,  in  matrix  notation. 


Z'®,.  0 


<E^(r, )E^tr_ ) >= 


11 

0 


e 

13 


0  0 
22 


\ 


e  0 

31  33 


/ 


From  (19),  an  identical  symmetry  holds  for  F. 


(20) 


With  this  information,  a  routine  calculation  following  that  described 


A  A 


in  [14]  shows  that,  for  b=  h  (see  (16)),  (18)  reduces  to 


,0.0  ^ 

Y^h  “  ©/(ttcos  9^)  dz' I  dz"A^^(z' ,  0)A(f,^(  z",  9) 


r.2.  ,  r  ff.  2  f  tf 

(sini<iF  (z  ,z  )+cos()iF  (z  ,z  )1 
11 


A 

while  for  b  »  v 


(21) 


Y^  =  It  (TTCOS  0  ) 
av  o  o 


dz' 


dz' 


App(z',9)App(z",e)[cos^<^  F^^(z'  ,z")+sin^<pF^^  (  z'  ,z")] 


(22) 


+cos()>[App  (z',9)A*  (z”,  8)F,,(z'  ,z")+A  (z',9)  A*  (z",0)F  (z'  ,z")] 

pz  13  pz  PP  31 

+  A  (z-  ,  e)A*  (z”,0)F  (z'  ,  z”)  > 
pz  pz  33  I 
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The  dependence  of  all  quantities  in  the  integrands  of  (21)  and  (22)  has 

been  suppressed  for  simplicity.  Eqs  (21)  and  (22),  in  conjunction  with  the 
B 

known  ,  provide  a  complete  description  of  the  scattering  coefficients 

(17).  It  is  interesting  to  observe  that  (21)  and  (22)  yield  non-vanishing 

cross  polarized  returns  in  the  special  case  of  backscattering  9=8  4»=’t. 

o, 

Of  course,  the  distorted  Born  approximation  does  not  lead  to  cross-polarized 
backscatter  cross  sections. 

An  approximation  which  is  very  useful  in  simplifying  later  calcula¬ 
tions  will  be  introduced  here.  For  a  sharply  peaked  correlation  function  c' 
in  (19),  the  expected  value  may  be  removed  from  the  integrand  to  yield 

F  (k,k  z',z")aE  ( z' )  w' ( z' ,  z",k  ,k  )  (23) 

kJl  X  y'  k£  x  y 

where,  from  (20), 

E,„(z')=  e  (r',r')  (24) 

k£  k£  ~ 

and 

w'(z',z",k  k  )=  dxdy  ( j x[ , j y | , j z-z ' | )exp[ ik  x+ik  y]  (25) 

X ,  y  /  X  y 

<.00 

Notice  that  translational  invariance  in  the  horizontal  direction  implies  that 
E  does  not  depend  on  x  '  or  y ' . 

Jv  X 

Equations  (21) -(25)  define  the  incoherent  wave  contribution  to  the 
scattering  cross  sections.  The  quantities  A^^,  App ,  etc.  have  been  dis¬ 
cussed  in  [5]  and  [14].  However  E  (z')  must  still  be  calculated  z'<0. 

k£ 

This  will  be  discussed  in  the  next  sedtion. 
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III  INTEGRAL  EQUATIONS  FOR  THE  INCOHERENT  FIELD  SECOND  MOMENT 

The  second  moment  dyad  is  governed  by  (11)  where  the  driving  term  is 

m  m* 

the  (assumed)  known  quantity  E  (r^)E  (r^)  exhibited  in  (13).  This  equation 

will  be  reduced  to  a  form  which  is  more  suitable  for  computational  work. 

The  principal  tool  in  this  reduction  is  the  representation  and  properties  of 
the  Green's  dyad  discussed  in  Appendix  A. 


When  eq.  (A-1)  is  substituted  into  (11),  several  delta  functions  arise 
which  can  iiiunediately  be  used  to  reduce  the  number  of  integrals  to  be 
performed.  It  is  found  that 


<E’^(r  )E'^*(r  )>-  (27r)"^ 
~  ~1  -  -2 


dz  dz” 


d^dh a'( 2 . z',  5  ,h)P(5,n,z',z")A' (z. 

ft  1  ®  .  ft  A 


•z",C,n) 


exp[iC(x^-X2)+in(y^-y2 )  ] 


(26) 


where  F  is  defined  in  (19).  If  the  correlation  function  c'  is  assumed  to  be 
sharply  peaked,  then  eqs.  (23) -(25)  will  be  applicable  for  use  in  com¬ 
puting  But  this  implies  that  the  solution  of  (26)  will  be  of  interest 

ab 

only  for  r  ■  r  -  r.  Thus 
"1  '*'2  " 


E(z)-  (2-rr)  “^[dz' dz"dCdnw' 


where  the  components  of  E  are  defined  by  (24)  and,  according  to  (13), 


<EE*>(  z  ')» 


m  ,  m  • 

E  (r')E  (r 
^  ^ 


')+E( 

ft 


z 


') 


(28) 
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Notice  that  the  first  term  on  the  right  hand  side  of  (28)  is  actually  a 

»  m 

function  of  z  only  because  E  depends  on  the  horizontal  coordinates  only 

through  a  factor  of  the  form  exp(  ik  sin9x*](see  [14]). 

o  o 

A  further  reduction  of  (27)  arises  from  the  observation  that  eq  (A-2) 
shows  that  the  angular  dependence  of  A*  occurs  as  a  simple  factors  of 
powers  of  cos<^  and  sin 41  if  cylindrical  coordinates  are  introduced  by 

means  of  the  equations 


t  *  pCOS(}l 

n*  0sin4)  (29) 

Since  the  correlation  function  was  assumed  to  be  spherically  symmetric, 
(25)  implies  that  the  ^,n  dependence  of  w'  (z' is  simply  a  depen¬ 
dence  on  the  parameter  p  only.  Performing  the  4>  integration  in  (28), 
recalling  the  symmetry  for  E  given  in  (20),  and  using  (A-2)  results  in  the 
equation 


11 


fO  f“  ( 

dz'  dz’  I  pdpW  '  <  Q 

»  *0  * 

00 


} 


E . Jdz'  dz"pdp W  +  Qjj<EE*  >^2  +  Qj^3<EE*  | 


22 


33 


(30) 


f°  f"  C  *  ) 

dz'dz"  pdpw'<Q^^<EE  ■>ii+23i<EE*>22+Q33<EE*>33[. 

ia,  •o  ^  ‘ 

Tdz'dz"!  pdpW 'S  P  <EE*>,,+  P  <EE*>,  ,  r 

j.  lo  <1  “2  ' 


E  -  E 
31  13 
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where 


Q^^(2,z'  ,2")=  (4n)  ^[a<^^(2,2'  )a^*^(2,2")+a^^(2,2'  )a^*<j  {  2 , 2",' 


-Q 


12 


Q^^(2,2  ,2  )=  (16Tr)  app{2,2  )-a(j,^  (2,z' )  app  (z,z")-a44,(z,2") 


Q^^(2,2' ,z")=  (4TT-)  ^ap2  (2,z' )a^2  {2,2") 


Q^^(2,2',z")=  (4ir)  a^p(2,2')a;p(2,2'') 


Q  (2,2',z")=  (2Tr)“^a,,(2,z' )a*  (2,2") 
33  ZZ  22 


P^(2,2',2")=  <4Tr)  |a^^(2,2' )+app  (2,z' )|a*^(z,2") 


(31) 


p^(z,z',z'0*  (4tt)  ^ap2  (z,z' )ap2  (z,z") 


The  dependence  of  all  quantities  in  (31)  on  the  parameter  P  has  been 
suppressed  in  the  notation. 


At  this  point,  it  is  useful  to  introduce  a  variable  change  in  the 
integrals  occurring  in  (30).  Notice  that,  according  to  our  assumptions, 
w'  decreases  very  rapidly  as  a  function  of  the  magnitude  of  z'-2"and  is 
nearly  zero  over  a  distance  beyond  a  few  correlation  lengths.  On  the  other 
hand,  the  z'  and  z"  dependence  of  the  remaining  terms  in  the  integrands  do  not 
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change  much  over  a  correlation  length.  Thus,  the  change  of  variables 
( compare  with  [ 5 ] ) 


$ 

z  =  z 
1 

u  =  z''-z'  (32) 


will  be  made.  With  the  above  assumptions 
argument  dependence  of  w'  as  w'(|u|;z^,p) 
that,  over  distances  of  a  few  correlation 


concerning  w*,  we  indicate  the 
where  the  z^  dependence  is  slow  so 
lengths  (|u|  not  too  large). 


W  (  u 


(Z^  +u,  p)»W 


!Z^,P) 


(33) 


* 

Liicewise,  <EE  >  does  not  change  much  over  a  correlation  length  so  that  it  may 
be  considered  to  be  a  function  of  alone  in  the  integrands.  Following  the 
procedure  in  [5]  and  using  these  assumptions  allows  (30)  to  be  approxi¬ 
mated  as 


fi  C  *  *  *  ) 

E  -I  dz^D,  ,<EE  >,,+  D,  ^<EE  b  <EE  > 

11  j  1|  11  11  12  22  Pz  33j 

E  -  dz  Jd  <EE*>  +  D  <EE*>  +  b  <EE*>  I 

22  1|  12  11  11  22  PZ  33 J 

E  =  dz  Jb  <EE*>  +  b  <EE*>  +2b  <EE* >  I 

33  1|  zp  11  zp  22  zz  33| 


(34) 


E,  dz  /d,  ,<EE*>,,+  l/2  C  <EE*>,*1 
13  I  11  13  13  PZZP  13 1 
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where 

3/4[bpp ]+  l/4bp(j, 


l/4[bpp +b^^-bp^  1 


(35) 


Sjpzz  ‘*’^4422  ^ 


Here 


°a6 


(z,z  )  = 
1 


-I 

(2Tr)  Re 


H 


i  ' 

duw 


(|u 


(36) 


where  the  subscripts  (aB)  are  one  of  the  pairs  {  pp  )  ,  (ijxj)),  (zz),  (pz)  or 
{tp)  and  Re  means  "real  part".  The  first  argument,  z,  in  the  components  of  the 
Green's  dyad  is  suppressed  for  simplicity  .  In  (35), 


bp^(2,Zi)=  (2Tr)  pdpjdu  w' [spp  ( z^)  ( z^ +u) +a,j„j,  (  z^  )a^  ( +u)]  (37) 


Also 


and 


1  r  * 

C^2ZP<^'2i)*  (2Tr)  |p<Spj^du  w'[a^^(z^)a^p(z^+u)+ap^(z^)a*p(z^+u)]  (36) 

(2t)-l^pdp  |du  w'[a^^(2^)a*^(z^+u)+a^„(z^+u)a*^(z^)]  '(39) 


where  the  subscript  o  is  either  P  or  4*  • 

m  m» 

Equations  (34  -(39)  are  a  set  of  integral  equations  for  E  when  E  g, 
is  specified.  However,  they  are  not  quite  satisfactory  for  numerical  studies 
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for  two  reasons.  First,  according  to  Appendix  A,  the  terms  in  the  kernels 
containing  a^^  have  an  explicit  delta  function  part  which  should  be  treated 
separately  in  numerical  work.  Second,  the  asymptotic  expressions  for  large 
values  of  the  parameter  P  given  in  Appendix  A  show  that  special  steps  must  be 
taken  in  order  to  avoid  numerical  difficulties  in  evaluating  the  integrals  in 
{36)-(39) . 

The  explicit  delta  function  contribution  from  terms  containing  a  is 

z  z 

easy  to  treat  and  results  in  an  explicit  evaluation  of  the  integrals  over 
in  eqs  (34).  Details  are  given  in  Appendix  B.  The  second  problem  concerning 
large  values  of  the  parameter  p is  treated  by  writing  all  integrals  over  p  as 


where  the  limit  P^  is  chosen  large  enough  so  that  the  components  of  the 
Green's  dyad  show  an  exponential  decrease.  A  reasonable  choice  is 

P^  *  k ^Max(ReJ^)  (41) 

where  the  maximum  is  over  the  range  of  z^  where  the  medium  is  random.  The 
asymptotic  expressions  developed  in  Appendix  A  may  be  used  as  an  approxi¬ 
mation  in  the  second  integral  in  (40).  This  allows  an  explicit  evaluation  of 
the  integrals  over  z^^  in  (34)  for  large p  when  the  order  of  integrations  is 
reversed  so  that  the  integral  over  (P^  ,»)  is  last.  The  details  are  found  in 
Appendix  C. 
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Of  course,  the  solution  of  eqs  (34)  requires  the  specification  of  the 
correlation  function  c'or,  equivalently,  its  transform  w\  We  will  write 
equations  under  the  assumption  that  C*  is  an  exponential  function  (it  will  be 
obvious  how  to  modify  the  following  results  for  a  single  exponential  if  a  sum 
of  exponentials  is  assumed  instead): 

2 

C'(lr-rj)=  <ie|  >exp[-|r-r,j/«]  (42) 

2 

where  <|^|  >  is  defined  by  (9)  and  (12)  for  r'=  r".  Then  W  (|u|;2)  is 
^  2 

W'«  <1^1  >V(|u|;p)  (43) 


where 


V=  {2'f/l)p  ^[l+p| u  |Jexp[-p|  u  iJ 


(44) 


and 


(45) 


The  slow  dependence  of  w'  indicated  by  its  second  argument  allows  the  possi- 
2 

bility  that  <|4|  >  and/or  i  may  be  a  function  of  z.  With  these  equations, 
all  remaining  integrals  over  P  in  the  interval  (p^,oo)  may  be  evaluated  com¬ 
pletely  as  is  shown  in  Appendix  D. 
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Upon  using  the  results  in  Appendices  B,  C,  and  D  in  eqs  (34),  it  is 
found  that  the  first  three  equations  of  the  set  (34)  may  be  written  as 


m  <EE 
2  22 


=  22*  "i^Ee\3  .£^jM33<EE* 


^  1  ^  <EE*>, 

11  11  22  p2  33 


E  +  m 
33  3 


<EE*>,,=  f  <3z,  <EE*>  +  f  <EE*>  +  M  <EE*>  t 

33  4  1  <  2p  11  zp  22  33  33  ) 


where 


M^^(z,z^).  3/4^^  +3/4f^^ +1/4^^ 

M^2(2,Z^)-  l/4[fp,+^^-^J 


(46) 


(47) 


M  (  Z ,  Z  )  :■  2f 
33  1  zz 

Here  the  subscripted  functions  f  are  identical  to  the  subscripted  functions  b 

defined  in  (36)  and  (37)  except  for  the  fact  that  the  upper  limit  in  the 

integration  over P in  (36)  and  (37)  is  replaced  by  P  (see  eq  (41))  and  a  is 

1  z  z 

cl 

replaced  by  a  (see  Appendix  B)  wherever  it  occurs.  The  multiplying 
factors  and  are  given  by 

■"l*  <141  >j-3P2  )/[32Jc^'‘|K^(z)|^]-3i'*P^(ip^)/32 

+  ll^F^(  Jlp^)Re[l/K^(z)]/[l6k^^  ]  I 

V  ^  *o‘  I  ^Pl  )/32 

(ip  )Re[l/K  (z)]/[l6k  ^]  [  (48) 
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"3=  < 


2  4  2,  rR  t  r° 


u  ;  2 ) 


(  z,2+u)+a'^^  (  z  ,z-ul/[jrk  ^  K  *(z)l  ( 
zz  o  o  ■*  i 


Here  the  functions  F  ,  F 

1  2 


and  F 
3 


are  defined  in  Appendix  D. 


All  of  the  functions  in  (46) -(48)  are  smooth  enough  so  that  no  numerical 
difficulties  should  result  in  the  solution  of  (46). 


In  the  same  way,  the  last  equation  in  the  set  (34)  becomes 


E^^+n<EES^^-JdzjN^<EE*>^^+  1/2  <EE*>*^  (49) 

where 


^13“  ^/^C%pzz  '''^4)$2z] 


(50) 


and  the  subscripted  g's  are  identical  to  the  c's  in  (38)  and  (39)  except 

that  the  upper  limit  in  the  P  integration  is  replaced  by  p  and  a  is 

1  zz 

cl 

replaced  by  a  .  The  multiplying  factor  n  in  (49)  is 
zz 


n-  7<(||2>/(8Jc^2)|-F2  ()lp^)./[)c^^|  K^(z)  ^|]+ f^F^  ()lp^  )/K^*  (z)  |  (51) 

+  f4'tr)c  (z)j  j  Pdp  du  W'{  |  u|  ;z)j  2[app  (z,z+u)+a^  ( z,  z+u)] 


+app  (z,z-u)+a^(z,z-u) 
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The  integral  equations  (46)  and  (49)  determine  the  second  moment  dyad 
E.  It  is  interesting  to  observe  that  (49)  is  not  coupled  to  (46)  so  that  the 
off-diagonal  components  of  E  may  be  determined  independently  of  its  diagonal 

elements.  Also,  notice  that  if  the  incident  wave  polarization  is  h^  (see  eq 

m  '' 

(31),  then  E  in  (28),  will  have  only  an  h  component  [14].  Thus,  the 
~  o 

driving  term  in  (49)  will  vanish  and  the  off-diagonal  elements  of  E  will  be 
zero  when  the  incident  polarization  is  horizontal.  On  the  other  hand,  all 
three  diagonal  elements  of  E  will,  in  general,  be  non-zero  when  (46)  is 
solved  for  either  horizontally  or  vertically  polarized  incident  waves. 


IV.  CONCLUSIONS 


A  system  of  integral  equations  in  one  independent  variable  has  been 
derived  for  the  second  moment  dyad  of  the  incoherent  field.  This  system  is 
well  suited  for  numerical  wor)<  because  all  singular  behavior  has  been 
extracted  from  the  integrands.  The  ]cernels  in  these  equations  are,  at 
present,  expressed  as  double  integrals  (see  (47)  and  (50)  together  with  (36- 
(39)).  However,  a  preliminary  examination  of  parts  of  these  Jcernels  has 
shown  that  one  of  the  integrations ( the  integral  over  u)  may  be  done 
analytically  if  the  techniques  discussed  in  [14]  are  followed.  Thus,  the 
numerical  wor)c  will  be  simplified  even  further  provided  that  the  remaining 
parts  of  the  Jcernels  can  be  treated  similarly. 

The  theory  developed  here  has  shown  that  a  calculationally  practical 
scheme  for  extending  strong  fluctuation  theory  beyond  the  distorted  Born 
approximation  for  scattering  problems  is  feasible. 
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APPENDIX  A.  GREEN'S  FUNCTIONS 


In  Section  II,  it  was  shown  that  the  Green's  dyad  F  for  the  operator 
(7)  was  required  in  order  to  determine  the  second  moment  of  the  incoherent 
field  within  the  random  medium.  As  in  [14],  it  proves  to  be  useful  to 
represent  F  as  a  two  dimensional  Fourier  integral.  The  discussion  in  [14] 
shows  that 


F(r,r')=  (21T  ) 


-2 


dC  dr)  A"  (  2,z' ,4,ri)exp|  i[4(x-x"  )-»n  ly-y if 
-L  '■  1 


(A-l) 


where 

(2  2 

<5/p'*zp 

2  2  2 

Here  p  =  %  +  t]  and  *pp »  functions  of  z,  z'  ,  and  p  .  Sub¬ 

stitution  of  (Al)  and  (A2)  into  L_r  =  6{r-r')  shows,  after  some  manipu- 
lation,  that  the  following  ordinary  differential  equations  are  satisfied  by 
the  components  of  the  Green's  dyad  in  the  region z<0; 


(|Ti/p^(app-a4x}>>  (|/p)3pz 


2  2 

a^^+(n  /p  )(app-a^^)  (n/plap^ 


(A-2) 


(q  /p )  a 


zp 


zz 


2 

d  ^ 

— 2  *  *  =  -  6(z-z  ) 


2 

^  ^PP 

dz^ 


b  ^!PP  +  a  app 


[a/(k, 


K  ) 
o 


J6(  z-z' ) 


(A-3) 


(A-4)  I 

I 
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(A-S) 

2 

^  ^pz  -  b^^pz  +  a  a  =  -  [ip/(k  )]  ^  ^  ^  +  [  ib/p]  5(  z-z ') 

di^  ^  o  o  Hz 


da. . 

a  =  [ip/a] _ Ph 

zp  dz 


(A-6) 


®zz  =  -  5(z-z-)/a 

where 

a(z)  =  k  (z)-  p^ 
o  o 


(A-7) 


(A-8) 


b(z) 


2  dK„ 

p  - /(aK  ) 

dz  ° 


(A-9) 


Continuity  of  the  tangential  components  of  the  electric  and  magnetic  fields 
at  each  point  of  discontinuity  of  (say,  at  z=-d)  requires  that  the  com¬ 
ponents  of  the  Green's  dyad  satisfy  the  conditions. 


j-d+o 

*4>4> 

'  -d-o 


=0 


(A-lOa) 


day“‘^'^°  =0  (A-lOb) 

dz “d-o 


‘PP 


-d+o 


-d-o 


=0 


(K  /.) 

o 


dz 


-d+o 


-d-o 


(A-lla) 


(A-llb) 
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I  -d+o 

ap  =0  (Al2a) 

I  -d-o 

j-d+o 

(K„/a)  _ Qz  =0  (A-12b) 

°  dz  -d-o 

These  continuity  conditions  also  hold  at z =0  where  K^(+0)=1.  These 
represent  the  boundary  conditions  at  z=  O.  There  is  also  a  condition  that 
the  components  of  the  Green ' s  dyad  behave  as  outgoing  waves  as  z  — > "®. 


It  was  shown  in  (31  that  it  is  very  convenient  to  work  with  the 

Riccati  equation  equivalents  of  (A-3)-(A-7)  if  the  additional  assumption 

(which  is  not  a  practical  restriction)  is  made  that  K  approaches  a  con- 

o 

stant  as  z — 9- a.  m  this  case,  one  may  define  functions  z )  by  the 
equations 


dz 


{A-13) 


where  a  ^  is  used  for  z>z'  and  a _  is  closed  for  z<z'.  The  boundary  condi¬ 
tion  atz'=0and  the  outgoing  wave  condition  as  z — imply 


0^(0)  ■  iV 


Q  (-«)  •  -iVa(-®) 

Substituting  (A-13)  into  (A-3)  (for  z^z’)  shows  that 


(A-14) 

(A-15) 


d  Cf+  2  _ 

_ ±  +  o±  +  a  =  0 

dz 


(A-16) 
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where  (from  A- 10b) 


(A-17) 


-d-o 


at  a  point  of  discontinuity  of  the  dielectric  constant  at  z=  -d.  Inte¬ 
grating  (A-13)  and  noting  the  delta  function  in  (A-3)  shows  that 


)|^|  z"  )dz''J/|o_(z ')-«^(2 ')J  i 

a_(z")dz''j4ja_{2  i 


if  2>Z 


(A-18) 


if  z<z 


For  some  purposes,  it  is  convenient  to  transform  (A-18)  by  means  of  the 

^  T 

reciprocity  relation  r(££')  =  r  (r,  r)(where  means  transpose).  This 
implies  that  an  alternate  expression  for  a^^is 


?|^-j  a_(z")dz'' j /|^a_(2)-a^  (z)J 


if  z>z'’ 


if  z<z' 


(A-19) 


in  a  similar  fashion,  functions  ^±{z)  may  be  introduced  by  the 


equations 


(3  . 

- »  P±  a„ 


(A-20) 
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These  satisfy 


A&i  +  P±-bp.  +  a  =  0 
dz 

where 

P+(0)  =  ia(0)/^K^(Q)>y/  j 


(A-21) 


(A-22) 


P_(-co)  =  -i<va(-®) 


{A-23) 


The  continuity  conditions  for  app  lead  to  the  conditions 


(K/a)Pdr  =0 
°  -d-o 


{A-24) 


for  Pi  at  a  point  of  discontinuity  of  the  dielectric  constant,  in  terms  of 
Pi  ,  app  is  given  by 


i(z  ")[^k^^Ko(2 ')j~^  exp[|  p+(z--)dz"]/Jp_(z ')-P+(z ')j  i 

z » 

i{z  ’  ):j^k^^K^(z' )j  ^  exp^  P -(  z”  ) dz"j/^P_(  z  ’  ) -P+(  z  '  )J  : 


if  z>z 


(A25) 


if  z<z' 


As  in  the  case  for  the  reciprocity  relation  leads  to  the  alternate 


representation 


(A-26) 


if  z<z' 
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The  component  ap^  of  the  Green's  dyad  satisfies  the  same  equation  as 
app  for  zj^z'.  Hence,  the  functions  may  be  used  to  describe  app.  How¬ 
ever,  the  delta  function  and  its  derivative  in  (A-5)  result  in  the  equation 


Thus, 


ap2  is  discontinuous  at  z=z^ . 


if  z>z' 


if  z<z' 


(A-27) 


Since  app  and  ap^  are  described  by  the  functions  ,  equations 

(A-6)and  (A-7)show  that  a„.  and  a  are  also  described  by  6+ .  Thus 

zp  zz  ^ 


^ipP_j_  (  z)/a(  z)  j  a 

^  (z)/a(2)j  a 


if  z>z' 


if  z<z'‘ 


(A-28) 


Like  Spjj,  this  is  also  discontinuous  at  z=z'.  Finally 


zz 


(z)/a(z)j  ap*  if  z>z' 

^ipP_  (z)/a(z)J  apz  if  z<z' 


(A-29) 


For  problems  where  integrals  arise  over  an  interval  containing  the  point 

z*z',  an  additional  delta  function  contribution  must  be  added  to  (A-29). 

This  arises  from  the  explicit  delta  function  shown  in  (A-7)  as  well  as  the 

da-_ 

fact,  evident  from  (A-5),  that  r*  has  a  delta  function  part.  The  total 

dz 

term  to  be  added  to  (A-29)  is  -ri/(k^  K  )1  6  (z-z'). 
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It  is  of  some  interest  to  note  the  asymptotic  behavior  of  the  com¬ 
ponents  of  the  Green's  dyad  as  the  parameter  p .  A  knowledge  of  this 
behavior  is  required  in  connection  with  integrations  over  intervals  con¬ 
taining  the  point  z=z'.  We  are  particularly  concerned  with  the  behavior  of 
terms  which  will  lead  to  non-zero  results  for  the  integrals  discussed  in 
Section  III.  Notice  that  (A-S)  and  (A-9)  show  that 

a(z)  {A-30) 

dK 

b(z)  -  — 2  /  K  (A-31) 

dz  ° 

Considering  first,  it  is  seen  that  the  boundary  conditions  (A-14)  and 
A-15)  reduce  to 

0^(0)  -  p  (A-32) 

0_(-oo)-*p  (A-33) 

Thus  (A-16)  yields 

(z)-*  +  p  (A-34) 
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independent  of  z  so  that  {A-14)  yields 

[l/(2p)]3xp  [-plz-z'U  (A-35) 

For  (A-22)  yields 

P(-«)-P  (A-36) 

and  (A-21)  yields 

PJz)-*p  (A-37) 

On  the  other  hand,  (a-22)  does  not  lead  to  P^(0)  -  p .  However,  as  a 

function  of  z,  the  solution  of  (A-21)  rapidly  approaches 

P+(2)  -  -  P 

with  the  thiclcness  of  the  transition  zone  approaching  zero  as  p 
these  results,  (A-29)  yields 

app  -►  -(  p/{2k^  K^(z') ) ]  exp  l-p  |  z-z' | ]  (A-39) 

These  results  immediately  imply 

apz  —  -(i  p  sgn(z-z')/(2k  (z'))]  exp  [-p(z-z'|]  (A-40) 

r  O  O 

and 

a  -•  -[  i  p  sgn(  z-z')/(2k  (z''))]  exp  [-plz-z'l]  (A-41) 

zp  O  O  ' 

2-28 


(A-38) 
-►  ®.  With 
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where  sgn(z-z')  is  the  sign  of  z-z'.  since  the  asymptotic  forms  of  a 


pz 


and  a^p  simply  change  sign  at  z=z',  integration  in  a  small  symmetrically 
placed  interval  containing  this  point  will  lead  to  zero.  Of  course,  if 
I z-2^|  is  not  small,  large  values  of  p  imply  that  integrals  containing  a 


pz 


and  a  will  approach  zero  because  of  the  exponential  factors.  Finally, 
zp 


zz 


[p/{2  K^(z'))]  exp  [-p  |z-z"|]  -  6(z-z^)/(k^  K^(z')) 


(A-42) 
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APPENDIX  B  ANALYSIS  OP  EXPLICIT  DELTA 
FUNCTION  PARTS  IN  KERNEL 


According  to  Appendix  A,  those  parts  of  the  integrals  of  eqs. 

(36)  and  (38)  involving  a^^  contain  an  explicit  term  with  a  delta  function. 
In  fact.  Appendix  A  gives 


cl  2 

a  (z,z,  )=  a  -6(z-z,  )/t)c  k  ) 
zz  1  zz  ^  o  o 


(B-1) 


cl 

where  a  is  a  function  in  the  classical  sense  containing  no  delta  functions 
zz 

and  is  given  by  eq.  (A-29). 


Using  (B-1)  allows  b^^(  eq  (36))  to  be  written  as 
pdpj  du  W(  lul;Zj^)  I  a°^(2^ 


b  =  (2Tr)~^  Re 
zz 


cl  ,  ,  cl*  , 

' “  )a  ( z  +u) 
zz  1 


■a^\z,  )6(z-z, -u)/(Jc  *  ( z )  I-a'^^*(  z, +u)  6(  z-z,  ) /[  )c  (z)J 

ZZl  O  O  zzi  loo 


+  6(  z-z^)6  (  z-zj. -u)/[lt^^| 


(B-2) 


,cl  -1 

ab  ~(2n)  Re 
zz 


pdp  /  W' ( I  z-z  l;z)a*^^(  z,  )/()t  *(z)]+  du  w'(lu|;z) 

,  (  1  zz  ^  CO 


a^*(z+u)6(z-z.)/(Jc  I-6(z-z,  )w' (0;z)/[)c  ^|K  1^] 

zz  *  o  o  ‘  o  1  o  I 

cl  cl 

where  b  is  given  by  (36)  with  a  replaced  by  a  and,  in  the  second  equa- 

Z  2  Z  Z  Z  Z 

tion  of  (B-2),  the  approximation  W'(|z-zJ  ;zj )aW' ( | z-z^ | ;z)  is  used.  When 
(B-2)  is  used  in  the  expression  for  E^^ (see (34)  ,  an  intergration  over  Zj 
is  performed.  We  consider  the  contribution  to  this  integral  from  the  sepa¬ 
rated  delta  function  parts  of  a  .  Thus 
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1-2 

j-co 


=  TT  ^Re  pdpJ  dz,  W' (  1  z-z  I  ;z,  )a  ^^(z  )<££*>,, (z)/[k  *(z)] 

Jo  I  r  1  '  1  zz  1  33  1  o  o 

^  (B-3) 

'0  cl  2 

+  duW'{lu  ;z)a  (z+u)<EE*>  (z)/[k  k  (z)] 

'  '  zz  33  o  o 

Laa 

-  W  (0;z)<EE*>,,(z)/[k  ‘^IK  (z)l^'li 
33  •-  o  I  o  !  Jf 

In  the  first  term  of  (B-3),  the  argument  of  <EE*>^^  may  be  replaced  by  z 
because  of  the  sharply  peaked  behavior  of  W  .  Further,  except  for  z  within 
a  few  correlation  lengths  of  the  surface  z=0,  the  upper  limit  of  the  first 
integral  within  the  curly  brackets  may  be  replaced  by  »  and  the  variable  of 
integration  changed  to  u=  z-z^.  Then,  making  use  of  the  fact  that  <EE*>^2 


is  real. 


I  =  d  <EE*>,,(z) 
zz  33 


fO  rO  cl 

2  dz,b  <EE*>,,(z  )=  2  dz  b  <EE*>,,(z)  -  d  <EE*>  ,(z)  (B-5) 

zz  33  1  1  zz  33  1  zz  33 

-30  tm 


where  the  multiplying  factor  d  is  given  by 

zz 


(iTk^  )  Re  pdp<j  du  W  ( I  ul  ;z)[2a^  (z+u)  +  a*^  (z-u)]/K  *(z) 
^0  '•Loo 


-W'(0;z)/[k^2|j^^^^,j2j^ 


(B-6) 


A  similar  treatment  of  a  in  (39)  shows  that  the  first  term  in  the 

zz 

last  equation  of  the  set  (34)  may  be  written  as 


^dz^Di3<EE*>^3(z 


l-l\ 


°13  <EE*>^3(z^)  -  dp^<EE*>^3(z) 


(B-7) 


cl 

where  D  is  the  same  as  0, _  except  that  all  terms  containing  a  are 
13  13  zz 


replaced  by  a 

zz 


2-31 


Report  9316 


The  multiplying  factor  dp^  in  (B-7)  is 


^6  =  I  pdp 


du  W* ( I u ■ ; 2 ) 


^00 


I  2[app  (z+u)+a^|^(z+u)]+  (  z-u) +a^(  z-u)  | 


(B-8) 
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APPENDIX  C  INTEGRALS  FOR  LARGE  PARAMETER  P 


When  eqs  (35)  -  (39)  are  substituted  into  (34)  and  use  is  made  of 
(40),  it  is  seen  that  integrals  of  the  form 


dz. 


Pdp 


du= 


dz. 


pdp 


^-00  ^  JOj  -1^ 


du  + 


dz. 


Pdp 


du 


(C-1) 


arise.  The  asymptotic  forms  for  the  compohents  of  the  Green's  dyad  given 
in  Appendix  A  may  be  used  as  an  approximation  for  all  of  the  a's  shown  in 
(36)  -  (39)  if  P>(^.  When  this  is  done,  all  integrands  contain  a  factor  of 


the  form 


a(  z,2^+u)  oc  exp[  -  P  iz-z^-u|  ] 


(C-2) 


where  the  subscripts  on  the  a  have  been  omitted  and,  in  contrast  to  eqs 
(34)-(39),  the  first  argument,  z,  has  been  included  because  of  its  role  in 
the  following  analysis. 

Consider  the  first  term  in  (C-1)  where  z>z^.  Eq  (C-2)  and  the  results 
of  Appendix  A  show  that 


a(z,z^+u)=  a( z,z^)exp( pu) 


(C-3) 


provided  that  the  approximation  K  (z,+u)a:K  (z  )  is  used  whenever  necessary. 

o  1  o  1 

For  large  p  ,  a(z,z^)  is  a  sharply  peaked  function  of  I  z-z^^l  .  Thus,  more 
slowly  varying  functions  of  z^  may  be  removed  from  the  integrand  when  per¬ 
forming  integrals  over  z^^.  Then,  denoting  any  pair  of  a’s  that  occur  in  the 
integrands  by  a^  and  we  find 


dz. 


pdp  a^(z,z^)a2*  (z,Zj^)F(p,z^)«j^pdpF(p,z)|  ^  ^ '  *3^ 


(C-4) 
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where  P(p,2^)  depends  on  the  integral  over  u  of  W  exp(pu)  and  the  components 
of  <EE*> 

On  the  other  hand,  for  2<2^,  W'(|u|;z^)  has  a  sharp  peak  at  u=0 
while  exp[-piz-2^-u|  ]  has  a  sharp  peak  at  u=  z-z^.  Therefore,  over  most  of 

i 

the  range  (z,0)  where  z^/z,  we  are  dealing  with  the  product  of  two  sharply 
peaked  functions  whose  peaks  are  separated.  Hence,  the  contribution  of  the 
second  term  in  (C-1)  is  negligible  compared  to  the  first. 

Thus,  the  integrals  in  (34)  for  which  p>p^  may  be  approximated  as 


dz. 


/•O 


I  dp  du  W  (  1  u.|  ;z  )a  (2,z  )a  *(  z,z  +u)<EE*>(2  ) 

0  I  1112  1  1 

J-oo 

00  (t) 


(C-5) 


a:|pdp[  du  W  ( |  u  [ ;  z )  e  [  dz  a  (  z ,  z  )  a  *  {  z ,  z  ]  <EE*>  (  z ) 
■la  Lx,  ^ 


* 

In  particular,  using  the  asymptotic  resutls  in  Appendix  A,  we  find  the 
following  for  integrals  over  the  range  (-<», z)  ; 

dZj^a^{z,z^)a,},^  (  z ,  z^  )a:l/[8P^] 

dz^app  (z,z^)ap^  (z,z^)»  dz^a^^C  z,z^ )  a^^*  (  z ,  z^  )a:p/[8k  (z)|^] 

’  -00 
E 

dz^app  (  z,z^)a^^*  (  z,z^)  s  -p/[8k^‘*|  K^(z)  | 
z 

dz^app(z,z^)a^^*  (z,2^)k  -^/[8k^^K*(z)p] 
z 

(z,z^)  *  l/[8k^  K^(z)p] 


(C-6) 


dziap^(z,z^)ap*(z,z,)  « 


dz  a  z,z,  )a  *  (z,z,  )  «  0 
1  zp  1  zp  1. 

^00 


cl 


In  (C-6),  a  is  the  part  of  a  which  does  not  contain  an  explicit  delta 
zz  zz  ^ 

function  (see  Appendix  B). 


2-34 


Report  9316 


APPENDIX  D  EXPLICIT  INTEGRALS  FOR  EXPONENTIAL 
CORRELATION  FUNCTION 


In  this  appendix,  the  completion  of  the  evaluation  of  the  integrals 
shown  in  (C-5)  will  be  made  assuming  that  the  function  w'  is  given  by  (43) 
and  (44).  The  integral  over  u  in  (C-5)  is  elementary  and  results  in 
o  . 

j  du  W'  eP*^  =  [2Tr  <  j^l  ^>/je]  p"^|[p+pj~^ +p[p-tp]“^|  (D-1) 

Substituting  (D-1)  in  (C-5)  and  noting  that  (C-6)  will  introduce  various  powers 
of  p  according  to  which  a's  occur,  we  find  that  three  )cinds  of  integrals  arise: 


(ip ) 
IF2  ) 
iF^  (ip  ) 

where,  upon  using  the  substitution  x  =  jp  in  the  left  hand  side  of 
is  found  that 


(D-2) 


(D-2) 


it 


Although,  it  does  not  appear  possible  to  express  the  integrals  in  (D-3) 
exactly  in  terms  of  elementary  functions,  numerically  satisfactory  approxima¬ 
tions  are  not  difficult  to  obtain.  First,  we  consider  limiting  behavior.  For 
large  X,  all  of  the  roots  in  (D-3)  may  be  expanded  in  powers  of  1/x  .  This 


leads  to 
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F^(X)  = 

1/(4X^) 

-  1/(4X^)  + 

= 

3/(4X) 

-  5/(12X^)  + 

F2(X)  = 

3/(20X^) 

-  5/(28X^ 

small  X,  it 

is  found 

that 

X 

II 

1  -  2X  + 

f  ♦  ... 

Fjdt)  - 

1-2^3 

3  3* 

*  ^1“  t  ... 

P3(X)  - 

2/X  + 

3  £n  X  +  . . . 

Using  (D-4)  and  (D-5)  as  guides  for  developing  functional  forms,  a  numerical 
calculation  of  a  rational  Chelbychev  fit  yielded 

Fj^(X)  =  (.83448  -  .71644X  +  .21851  X^)  /(. 83473  +  .93764X  +  X^) 

F2(X)  =  f  -  (.87871  -  .58270X  +  .  15574X^ )/( 1 . 31823  +  .59901X  +  X^) 

F2(X)  =  2/X  +  3£nX  -  (3.1544  +  .6467X  +  5. 0162X^)  /  (2.9216  +  .6307X  +  X^ ) 

(D-6) 

for  X  <  1  while 

F^(X)  =  (.25003  X^  -  1.375  x  lo”^  )/[  X^ (X+^ ,94783  )] 

F2(X)  =  (.74972  X^  +  .13296)/  [x  (X^  +  .72405)]  (D-7) 

D^(X)  »  (.15010  X^  -  .01927)/  [x^(X^  +  1.07866)] 

for  X  >  1. 

Comparison  of  numerical  results  from  (D-6)  and  (D-7)  with  high  preci¬ 
sion  numerical  evaluation  of  (D-3)  showed  that  the  rational  approximations  to 
F^,  F^,  F^  have  maximum  errors  of  .06%,  .03%,  and  .06%,  respectively. 

The  integrals  in  (C-5)  are  completely  specified  by  F^^,  F^f  F^  and  the 
various  numerical  factors  shown  in  (D-1),  (D-2),  and  (C-6). 
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